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ABSTRACT 

The satellite COROT will search for close-in exo-planets around a few thousand stars 
using the transit search method. The COROT mission holds the promise of detecting 
numerous exo-planets. Together with radial velocity follow-up observations, the masses 
of the detected planets will be known. 

We have devised a method for predicting the expected planetary populations and 
compared it to the already known exo-planets. Our method works by looking at all 
hydrostatic envelope solutions of giant gas planets that could possibly exist in arbitrary 
planetary nebulae and comparing the relative abundance of different masses. We have 
completed the first such survey of hydrostatic equilibria in an orbital range covering 
periods of 1 to 50 days. 

Statistical analysis of the calculated envelopes suggests division into three classes of 
giant planets that are distinguished by orbital separation. We term them classes G 
(close-in), H, and J (large separation). Each class has distinct properties such as a 
typical mass range. 

Furthermore, the division between class H and J appears to mark important changes 
in the formation: For close-in planets (classes G and H) the concept of a critical core- 
mass is meaningless while it is important for class J. This result needs confirmation 
by future dynamical analysis. 

Key words: planets and satellites: formation - planetary systems: formation 



1 INTRODUCTION 



Since the discovery of 51-Pegb (|Mavor fc Quelozl Il995f ). 
more than 200 exo-planets have been discovered, most by 
the radial velo city technique. This year, the satellite mis- 
sion COROT (|COROTl [2006) will be launched hoping to 
add many more planets to the list. The COROT satellite 
will be on a two-fold mission: It wil l do A) astroseismology 
(|Baglin fc The COROT Teamlll998h and B) look for planets 
using the transit search method dBorucki fc Summersll 19841 ; 
ICharbonneau et all l2000l : iRouan et al l I2000T ). The transit 
search programme hopes to find a relatively large number of 
planets. It is the task of theoreticians to make a prediction 
beforehand. 

The standard giant planet for mation model is the so- 
called core-accretion model as in iMizunol (|l980h . In this 
model, planet formation starts with sedimentation and co- 
agulation of the condensible material i nto small solid cores 
jGoldreich et all 120041 : iLissaueij Il993l : IWetherill fc Steward 
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GUI). As this core grows, it becomes massive enough to 
gravitationally bind some gas. Consequently, it acquires an 
envelope of gas and dust. The evolu tion of this envelope has 
been studied by many au t hors (e. g . Bodcnhcimcr fc Pollack! 
1 19861 : iPollack et all 19961 : IWuchterl|[l990l . Il991al lbh. 

It should be mentioned, that planets could also form 
as de scribed by the gravitational instability scenario (|Bossl 
2002). Nevertheless, today's planets are in better agreem ent 
with the core-accretion scenario (see lSantos et al]|2005l ). In 
this paper we work on the basis of the core-accretion sce- 
nario. 

A natural procedure when trying to predict the distribu- 
tion of giant planets is the statistical approach: Calculate the 
evolution of a large number of randomly placed planetesimal 
" seeds" , starting with small planetesimals and letting them 
evolve to the final planet. For each seed, the full evolution 
is calculated including core growth, accumu l ation of enve- 
lope, m igration, etc. (see e.g. lBenz et al.|[2006l ) (lAlibert et alj 
l2005dlblldY This has the advantage, that a large number of 
processes can be included into the algorithm. Furthermore, 
the physical evolution is modeled in a natural way. However, 
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there are disadvantages as well. First of all, these calcula- 
tions are computationally intensive and a very large number 
of such calculations need to be done in order to gain a sta- 
tistically significant result. In addition, one needs to know 
the exact environmental conditions in which the evolution 
of the seeds should be calculated. This is a problem. Even 
in our own solar system, the conditions during the time of 
formation are only vaguely known. Nebula densities range 
between two extremes: There must have been enough mate- 
rial to form all the planets (the concept of the min i mum mass 
solar nebula: lHavashi|[l98ll : lHavashi et ailll977l . Il979l ^ and 
the nebula must be gravitationally stable. F or a more thor- 
ough discussion on the nebula variety see IWuchterl et al.l 
(2000). For other stars, the primordial proto-planetary neb- 
ula is constrained even less. As long as these important pa- 
rameters are not known with some degree of precision, we 
think that it will be difficult to make a good prediction in 
this way. 

Therefore, we use a different approach. We study all 
possible equilibrium states consisting of a solid core and a 
gaseous envelope with cores of different sizes and a range of 
nebula densities. 

It is our goal to make a prediction for COROT. As this 
satellite mission will only be sensitive t o planetary orbits 
shorter than 50 days l)Borde et al.ll20o"3 ). we'll restrict our 
prediction to close-in planets ranging from I to 50 days or- 
bital period. Our prediction is only valid for gas giant plan- 
ets, terrestrial planets' mass distributions cannot be pre- 
dicted in this way. 



2 PLANET PREDICTION METHOD 

Using a wide range of nebula pressures and core masses 
we can calculate the possibly existing envelope-core combi- 
nations in hydrostatic equilibrium. Assuming that all such 
states are equally likely, the relative frequency of planetary 
masses (core+envelope mass) corresponds to a distribution 
function of planet masses. 



2.1 Calculation of the envelope structures 

Before we discuss the prediction method in detail, we'll de- 
scribe how the individual planet-envelope structures are cal- 
culated. 

Each "planet candidate" consists of a solid core of 
fixed densitjQ that is embedded in a nebula of a par- 
ticular pressure. The envelope structure is determined 
by the well-known equatio ns of stellar structure (e.g. 
iKippenhahn fc Weigertl Tl990). These, we calculate in radial 
symmetry and neglect rotation. The effect of rotation i s neg- 
ligibl e in all but the most extremely rotating cases l|Gotj 
1 19891 ). 

A constant infall of planetesimals onto the core releases 
gravitational energy that is transported through the enve- 
lope by either radiation or convection. We use the diffusion 
approximation for the radiative energy transport and zero 
entropy convection. 

The properties of the envelope are determined by the 



we use the value of p = 5500 kg m 



equation of state bv lSaumon et all (|l995l ). Rosseland-mean 
opacities ft(p, T) are interpolated from a combined ta- 
ble: Opacities includ e Rosseland-me an dust opacities from 
Pollack et ail l| 19851 . lgT sC 2.3), Alexander fc Ferguson! 
( 1994 ) values in the molecular range, and lWeiss et al.1 (119901 ) 
Los Alamos high temperature opacities. 

The planet extends out to the hill-radius where the pres- 
sure of envelope an d nebula are set to be equ al. 

As discussed in iBroee fc Wuchterll |2006), the problem 
is fully specified when the following six quantities are spec- 
ified: 
The 

(i) core mass 

(ii) pressure at the core P C ore, 

(iii) mass of the host star M*, 

(iv) semi- major- axis of the planet a, 

(v) nebula temperature T nc b, and the 

(vi) planetesimal accretion rate M . 



The parameters (i) and (ii) are our independent parame- 
ters. By varying these two independent parameters, we can 
determine all possible hydrostatic envelope solutions for a 
given "location". A "location" is determined by the param- 
eters |(iii)|(vi)| They give the environmental conditions of 
the proto-planet. 

Parameters (iii) | (iv) and (v) are determined by the host 
star and the location of the planet. The nebula temperature 
can be calculated in thermal equilibrium with the star: 



Tneb — 280 - 



_a_\-i/2 / L, 
1AU/ \L Q 



1/4 



K 



(1) 



with L* the lumino sity of the planet host star and the so- 
lar luminosity (see|Hayashi 1 98ll : lHavashi et al.| [l985) . This 
implies a passive disk, i.e. no viscous heating and assumes 
that the nebula is optically thin. 

The only remaining free parameter is the planetesimal 
accretion rate M. Proper values in agreement with plan- 
etesimal theory range from M — 10~ 2 a" 1 (note that a 
stands for one year here, not the semi-major-axis) very close 
to the star to M = 10 -6 Mg^ a" 1 at Jupiter distance^. 

For a detailed description of the equatio ns and bound- 
ary conditions, see lBroeg fc Wuchterll ([2006 ). 



2.2 Calculating a mass spectrum for a fixed 
location 

2.2.1 A set of solutions for a fixed location: the manifold 

For a given location as defined by the parameters (iii)|(vi) 
in section 12.11 we can calculate all hydrostatic equilibrium 
solutions to the equations of stellar structure. One such 
set of solutions co vering a wide range in th e M c ,P c -plane 
we term, following IPecnik fc Wuchterj (|2005l ). a " manifold" 
or "solution manifold" for the given location. Each mani- 
fold contains, once calculated, all envelope structures that 
can possibly exist hydro-statically inside any nebula at the 
given location. One example for such a manifold is given in 

2 The high value corresponds to a orbital distance of a = 0.04 AU, 
particle-in-a-box planetesimal accretion theory with a gravita- 
tional enhancement factor _F g ?a 20, a minimum mass solar nebula 
[|Havashilll98ll; lHavashi et al.lll985l) and Jupiter mass objects. 
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Figure 1. Manifold of proto- planet masses for a four day orbit 
around a solar type host star. The accretion parameter is set to 
M = 10 — 4 a -1 . The total proto-planet mass as a function of 
the parameters M c , P c . All axes are logarithmic. The results are 
connected to show a 3-dimensional surface. This surface is colored 
using the color-map shown in the upper right corner mapping the 
outside pressure P x in a logarithmic way. This shows that the out- 
side pressure is varied from fa 10 4 fo fa 10~ 30 Pa. Everything with 
higher nebula pressures is by construction gravitationally unsta- 
ble, lower pressures correspond to vacuum and cannot be calcu- 
lated here properly because we do not calucate an atmosphere. 
This plot covers more than 6 magnitudes in M c and Mtot , and 5 
magnitudes in P c . 

figure [1] It shows the total mass of the proto-planet as a 
function of the parameters M c and P c . 



2.2.2 Deriving the mass spectrum, from a manifold 

Having calculated a manifold, we now make the following 
assumptions: 

(i) All equilibria are equally probable. 

(ii) All equilibria are stable and can be dynamically 
reached. 

Now we can - quite in analogy to statistical mechanics 
- derive a distribution function for various properties of the 
proto-planets at that "location". The quantity we are inter- 
ested in is the mass of the planet. By quite literally count- 
ing off the occuring masses in the manifold we can derive 
what we call the "mass spectrum": The relative frequency 
of planet masses at this location^. 

For this application of the manifold, it is important not 
to choose a certain range of core masses and core pressures 



3 In order to produce a histogram of continuous data, the data 
have to be binned to a fixed bin size. We chose a logarithmic 
binning with a bin size of 0.05 dex. 
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Figure 2. Mass spectrum of proto-planet masses for a four day 
orbit around a solar type host star. The accretion parameter is 
set to M = 10 -4 M ffi a -1 . The area is normalized to 1. The red 
lines mark the value of 1 Mq and 1 M-^ . The mass-distribution for 
IMq at 4 d is clearly dominated by two peaks, one at fa 17Mgj 
and another at fa 210MQ or fa 0.6 M^. 

implicitly by choice of a scale. Therefore, we chose a scale- 
free distribution to sample the parameter space in M cole 
and Pcorc- In this way no core mass or core pressure is se- 
lected and all values are treated alike. The only scale- free 
distribution is a power-law distribution, or correspondingly 
a log-equidistant sampling of the parameter space. Obvi- 
ously, if there were a dominant scale, e.g. in the core mass, 
this would drastically change the outcome as compared to 
our scale-free set-up. 

The result for the same location as figure [1] is shown in 
figure [2] This mass spectrum is derived based on assump- 
tions (i,ii). Therefore, either dynamical instability of some 
parts of the manifold, or quasistatic contraction with signifi- 
cant mass gain will change the observed mass spectrum. On 
the other hand, agreement of the mass spectra with obser- 
vation would be a strong indication for the validity of the 
assumptions (i,ii). 



3 PLANET PREDICTION RESULTS 
3.1 Manifold survey 

We have calculated manifolds and corresponding mass spec- 
tra for a wide range of locations by varying the following 
three parameters: 

(4)T orb : 1, 4, 16, 64 days 

(4)M» : 2, 1, 0.8, 0.4M Q ; L, = 16, 1, 0.42, O.O4L 

4 This corresponds roughly to spectral types A2, G2, Kl, or. M2. 
Luminositie s of the host star are assigned to the masses following 
iGravl Jl992l) . 
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(3)M: 10" 2 , 10~ 4 , lO^M^a -1 . 

This results in a 3-dimensional grid of locations, a total of 
48 manifolds. This is the first complete survey of hydrostatic 
proto- planets i n close orbits. The full set of results can be 
scon m iBroed (|2006b) and is also available on-line at Broeg 
(2006aj). This survey - named Corot survey Mark 1 vl.l - 
revealed a large diversity of mass spectra in the range from 1 
to 64 days orbital period. The host star mass also has large 
impact on the mass distributions. 



3.2 Statistical properties — Three classes of gas 
giants 

As stated in section [2.2.21 a manifold can be used to deter- 
mine the mass spectrum using the following two hypotheses: 
1) all equilibria are equally probable and 2) they can be dy- 
namically reached, i.e. there exists a track from some set of 
initial conditions to each state. Using these hypotheses we 
can derive several interesting properties of the giant planets. 

One major result of this survey is the fact that all mass 
spectra for close-in orbits exhibit two peaks. This is so for all 
tested values of M. Moving to larger orbital distances, these 
peaks move closer together and eventually merge into one 
peak. For a solar type host star, this happens at an orbital 
period of around 16 days. 

The full set of mass spectra of our survey leads to the 
grouping of the planets into three classes: 

Class G Extremely hot gas giant jf| reside very close to the 
host star. Their surface temperature is above dust sublima- 
tion temperature. Planets in this class have a very large 
upper mass limi10. For solar type host stars, the mass limit 
is roughly at 2.5M-2 f Q For host stars of 2Mq, this limit 
is extended up to 6M^. More massive planets should not 
exist this close to the star. In addition, we expect a very 
large quantity of so-called hot Neptunes with masses around 
16 Mq corresponding to a large second peak in the mass 
spectrum. 

Class H Hot gas giants reside in-between the classes G 
and J. Their surface temperature is below dust sublimation 
and they are close enough to their host star so that the mass 
spectra still show two distinct peaks. We expect them to be 
less massive than 1 M^ (for a 1 Mq host star). 

Class J Jupiter-like gas giants show only one peak in their 
mass spectrum. Class J planets can be much more massive 
than the classes G and H because the equilibria can gain 
significant amounts of mass by quasi-static contraction while 
the nebula is still present (see section [3~3j) . 

For a solar type star, the boundaries between the groups 
G,H and H,J are at 4 days and 16 days orbital period, re- 
spectively[f| 



5 German Ganz heifi 

6 derived as the largest occuring masses in the mass spectrum 

7 This value is strongly dependent on M. The given value corre- 
sponds to M = 10 _4 Mq a -1 . It is 5.3 and 0.8 for an M of 
10~ 2 and 10~ 6 respectively. 

8 The 16 day boundary depends strongly on the host star mass 
and the planetesimal accretion rate. 16 days correspond to a 1 Mq 
star and a low accretion rate (M = 10 — 6 Mq a -1 ). For a slightly 



3.3 Discussion 

As discussed in the last section, the transition from class 
H to J is marked by the merger of the two peaks in the 
mass spectrum. We have performed b oth isothermal linear 
instability analysis (Schonkc, J. 2005) and isothermal non- 
linear instability analysis (|Pecnik. B.II2005T ) of a number of 
manifolds. These calculations suggest a fundamental change 
in dynamical properties that coincides (or is caused by) the 
merger of the peaks: In the one-peak case, entire regions 
in the manifold appear to be unstable owing to transitions 
between two states of similar mass. 

Another change in behaviour at the merge position can 
be derived from the manifolds directly: At large orbital dis- 
tances, there is a well-defined critical core mass beyond 
which no more static solutions exist inside a nebula. The 
value of the critical core mass depends only weakly on neb- 
ula pressure. This leads to the explanation why Jupiter and 
Saturn appear to have similar cores. 

At small orbital distances, however, the critical core- 
mass becomes very strongly dependant on nebula pressure. 
In consequence, it is always possible to find a nebula pres- 
sure where the core is sub-critical, i.e. where a static solu- 
tion exists. This renders the concept of a critical core mass 
meaningless for classes H and G. Following the above line of 
reasoning for close-in planets, it follows that no significant 
mass gain is expected by the disappearance of the nebula. 
Therefore the calculated mass spectrum of classes G and H 
could be very much like the observed mass spectrum in this 
regime. 

These considerations have been tested in a small num- 
ber of calculations using full radiation hydrodynamic planet 
formation. One such case is the planet HD 149026 b at a dis- 
tance of 0.042 AU corresponding to an orbital period of 2.87 
days. Ou r calculations show a com pletely hydrostatic evolu- 
tion (see lBroeg fc W uchtcrl 2006) without significant mass 
gain in the final phase. The final planet has the same mass 
as the equilibrium configuration in the manifold. 

3.4 First comparison with observations 

ISantos et all (|2005l ) detect a paucity of high-mass planetary 
companions with orbital periods shorter than ~ 40 days. 
This is in agreement with our separation in upper mass limit 
classes (G&H) and the J class without a strict upper mass 
lirmL 

iGaudi et al.l (|2005l ) went a step further dividing the ob- 
served exo-planets into very hot (VHJ) and hot Jupiters 
(HJ) with a dividing line at 3 days orbital period. They 
observe that the VHJ exhibit higher masses than the HJ; 
Specifically, the VHJ masses are larger than lM^. This is 
in perfect agreement with our separation into groups G & 
H and the upper mass limit of 1 for group H. 

As a final comparison of our method to observations 
we performed a direct "prediction" for the host stars of to- 
day's exo-planets. Because of the class properties, nameley 
the stability of equilibrium states for classes G & H and the 
lack thereof for the J class, we only compare our predicted 



higher accretion rate (M = 10 4 Mq a 1 ), the two peaks in the 
mass spectrum merge around 32 days orbital period. 
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Figure 3. Predi cted and o bserved mass distributions (Observed 
exo-planets from Schneider 2006, on July 7, 2006). For this pre- 
diction we used only solar type host stars. The high-mass peak 
is in good agreement with the observed data. Only the very high 
data point is not reproduced but it might be a brown dwarf with a 
different formation mechanism. Including higher mass host stars 
would produce slightly higher upper masses in our prediction. 
The lower-mass peak of our calculations can already be noticed 
in the observed data but these planets are at the detection limits. 
We expect this part of the observed distribution to grow as the 
instruments become more sensitive. The high-mass end, on the 
other hand, should be complete for the observed stars. 



4 CONCLUSION 

We have presented a new method to predict the mass dis- 
tribution of gas giant planets that analyses all possible hy- 
drostatic equilibria. It has the advantage of not needing the 
nebula density as input. Only nebula temperature and plan- 
etesimal accretion rate must be known. In a passive disk, this 
can be easily approximated using the host star properties. 
This leaves the planetesimal accretion rate as the only free 
parameter. 

Using our new method, we are able to split the giant 
planets into three classes G, H, and J which have distinct 
properties (see section [33}. We compare these properties to 
the observed mass distribution of the exo-planets and find 
good agreement. We also produce a mass distribution for the 
exo-planet host stars having close-in planets and can repro- 
duce the observed mass distribution. The agreement with 
observations is a strong argument that the equilibria are in- 
deed dominating the formation process of close-in planets 
and that a large variety of proto-planetary nebulae is in ex- 
istence. 

To produce a prediction for COROT star fields, the ex- 
isting mass spectra have to be averaged according to the 
distribution of stars in the COROT fields and a concept 
has to be developed to determine the relative planet abun- 
dance of planets at different orbital distances. So far we can 
only predict planetary mass distributions at given orbital 
distances. 



masses at orbital periods less than 20 days. At the time 
of our analysis, 54 exo-planets fell in that regime. For the 
prediction we assumed all host stars to be solar(f] The ref- 
erence mass spectrum was obtained by binning the masses 
of these 54 exo-planets in 0.3 dex mass bins. As most of the 
detected exo-planets have been observed using the radial ve- 
locity technique, we added a 30% correction to the observed 
planet masses as is statistically expected. This we compared 
to our theoretical mass spectrum: For each detected exo- 
planet we computed a mass spectrum using the correspond- 
ing orbital distance. All such mass spectra were added and 
binned according to the reference mass spectrum0 The re- 
sulting mass spectrum is compared to the reference spec- 
trum in figure [3] Please note that our calculated mass spec- 
trum can only be expected to reproduce the observed one if 
there are no phases of quasi-static contraction while the neb- 
ula still exists, i.e. while there is a significant mass reservoir 
for the planet (this is the case for HD 149026b). 



9 Most exo-planet host stars to date are of solar metallicity. 

10 More precisely, we used the following procedure to obtain 
the theoretical mass spectrum: 1) The detected exo-planets are 
grouped into bins of orbital period (0—2,2-8,8-20 days). 2) The 
computed mass spectra for the planetesimal accretion rates M = 
10" 2 and 10" 4 M e a- 1 at a given orbital period (1, 4, 16 days) 
are multiplied with the number of exo-planets in the correspond- 
ing period range and added together. 3) The resulting mass spec- 
trum is then normalized and binned like the reference spectrum. 
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